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Abstract 

The structure and the degeneracy of the ground state of a fully-frustrated 
XY-model are investigated for the case of a dice lattice geometry. The re- 
sults are applicable for the description of Josephson junction arrays and thin 
superconducting wire networks in the external magnetic field providing half- 
integer number of flux quanta per plaquette. The mechanisms of disordering 
of vortex pattern in such systems are briefly discussed. 
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I. INTRODUCTION 



De Gennes |IJ and Alexander B have shown that the linearized Ginzburg-Landau equa- 
tions for a superconducting wire network in external magnetic field can be mapped on the 
eigenvalue equations for a single electron hopping problem in the same geometry. Thus 
the dependence of a (mean-field) superconducting transition temperature on external field 
can be found by following the field dependence of the lowest eigen-value in single electron 
problem. 

Recently it has been shown by Vidal et al |3j that the single electron hopping problem 
has very special features in the case of so-called dice lattice (Fig. 1) if the value of 
the magnetic flux per plaquette $ is equal to one half of the flux quantum <3>o- Namely, the 
spectrum of an electron lacks any dispersion and is reduced to three discrete levels. 

This theoretical result has led to the interest in experimental investigation of supercon- 
ducting networks with dice lattice geometry |||7]]. It has been shown [0] that in contrast 
to the case of $ = $o/3, for which vortex pattern in dice network is nicely ordered, at 
$ = $o/2 [so-called fully- frustrated (FF) case] the vortices do not form (presumably at the 
same temperature) any regular pattern. The authors of Ref. J7J have suggested that this ab- 
sence of ordering is related with an infinite degeneracy and localized structure of the states 
corresponding to the lowest energy level in terms of a single electron problem (or to the 
lowest free energy in terms of a superconducting network). 

Although this conjecture can be correct, it still has to be verified. A superconducting 
network problem reduces to a single electron hopping problem only at the mean-field tran- 
sition point. Below it the non- linear terms in Ginzburg-Landau equation become important 
and may completely or partially remove the high degeneracy of the state with the lowest 
free energy. 

In the present work we use a different approach for theoretical investigation of an or- 
dering in a FF superconducting system with dice lattice. We consider another limit when 
the amplitude of the order parameter is well defined and uniform, but phase fluctuations 
are possible and can lead to destruction of an ordered state. In that limit a discrete super- 
conducting system (a wire network or a junction array) in external magnetic field can be 
described by a frustrated AF-model introduced in Sec. 2. 

In Sec. 3 we propose a highly symmetric state which due to simplicity of its structure 
may be a good candidate for the ground state of a FF Xy-model with dice lattice. In Sec. 
4 we show that this state has high additional degeneracy because it allows for formation of 
zero-energy domain walls. Sec. 5 is devoted to a brief discussion of possible consequences 
of this additional degeneracy for the disordering of vortex pattern in FF superconducting 
systems with dice lattice geometry. 

II. THE MODEL 

In the regime when only phase fluctuations are of importance an array of weakly coupled 
superconducting islands can be described by the Hamiltonian: 

tf = £v(%), (i) 

(ij) 
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where the sum is performed over all pairs of coupled islands and 

Oij = <fj -Vi- 1 ^- f dx A(x) = -9ji (2) 

is the gauge-invariant phase difference which can be associated with the link (ij). Here ipj is 
the order parameter phase of j-th superconducting island and A is the vector potential. The 
phases tpj are defined up to a shift by a multiple of 2tt, therefore the interaction function 
V(9) has to be periodic in fl. The form of V{6) depends on the type of the coupling. For 
Josephson junction array 

V(6) = -Jcos6, (3) 

where J is the coupling constant of the junction. 

Summation of Eq. (||) over a perimeter of a lattice plaquette imposes a constraint 

E% = - 2 ^ ( 4 ) 

□ 

where the frustration parameter / is equal to $/$o an d $ is the magnetic flux threading 
the plaquette. In the limit when screening effects can be neglected $ is determined by the 
external field and for the uniform field and flat geometry is proportional to the area of the 
plaquette. If all the plaquettes have equal areas the value of / is the same for all plaquettes. 
In such case a system is called uniformly frustrated. 

When interaction function V(9) is periodic in fl it is convenient to consider variables fly- 
reduced to the interval [— ir, it]. That makes the description of any state in terms of fl^- more 
transparent, but transforms the constraint (||) into 

J20ij = 27rM, M = m-f; (5) 

□ 

where m is an integer. The form of Eq. ([]) shows that / can be reduced to the interval 
— 1/2 < / < 1/2, all other values of / being equivalent to some value from that interval. If 
V(fl) is an even function / is additionally equivalent to — /. The case of / = 1/2 is usually 
called a fully-frustrated (FF) XF-model. 

Well below mean-field transition temperature a network of thin superconducting wires 
can be described by the same Hamiltonian (P with the term ipj — ip^ in the definition of fl^- 
[Eq. (0)] now substituted by the integral // dx(dip/dx) along the link (ij) ||. In that case 
summation of Eq. (fj) around a perimeter of a plaquette leads directly to Eq. (H). In the 
hmit of long thin wires the interaction function is almost harmonic M : 

V(6) oc 6 2 . (6) 

In this work we use the term "frustrated Xy-model" for a system defined by Eqs. flU) 
and (|5|) with a general form of the interaction function and not only for V(9) of the form 
Thus our approach is valid for the description both of junction arrays and wire networks. 
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III. THE GROUND STATE 



In a FF XF-model all variables M are half-integer and different low-lying extrema of the 
Hamiltonian can be characterized by the distribution of positive and negative half-vortices 
(M = ±1/2) in the plaquettes of the lattice. The vortices of the same sign repeal each other, 
therefore in the ground state they can be expected to be situated as far from each other as 
possible. In particular, in the case of a FF XF-model with square lattice the positive and 
negative half- vortices form in the ground state a regular checkerboard pattern ||. Analogous 
pattern in which the nearest neighbors of each half- vortex are of the opposite sign is possible 
when they occupy the sites of a honeycomb lattice, that is in FF AF-model with triangular 



A dice lattice is dual to a Kagome lattice, therefore in a FF AF-model with dice lattice 
the half-vortices can be considered as occupying the sites of a Kagome lattice. Since a 
Kagome lattice is constructed from triangles it is impossible to distribute the half-vortices 
in it in such a way that all nearest neighbors are of the opposite sign. The half-vortices of 
the same sign will have to form clusters and the minimal size of such clusters which allow for 
covering of a Kagome lattice turns out to be equal to three. The most symmetric example 
of a regular (periodic) arrangement of vortices on a Kagome lattice in which half- vortices of 
the same sign form the clusters of the size three (triads) is shown in Fig. 2a. 

This state has the 12-fold degeneracy and can be described as a regular lattice of vacancies 
(absent positive vortices) on the background of / = 2/3 ground state or, equivalently, a 
regular lattice of extra positive vortices on the background of / = 1/3 ground state (cf. 
with Ref. 0]). However symmetry considerations show that its structure in terms of gauge- 
invariant phase variables 9ij (which is shown in Fig. 3a) is very simple and can be constructed 
by repetition (with rotation and reflection) of a simple three-link pattern shown in Fig. 3b. 
Although we can not rigorously prove that this state has the lowest possible energy (which 
is typical for AF-models with nontrivial frustration), we believe that the simplicity of its 
structure strongly supports this conjecture. 

In the state depicted in Fig. 3a the variables 9ij acquire only three different values which 
we denote 9 a (a = 1, 2, 3; < Q\ < 9 2 < 9 3 < it) and show in figures as single, double and 
triple arrows. The half-vortices of the same sign are separated by single arrows, the central 
half-vortex of each triad is separated from its neighbors of the opposite sign by triple arrows 
and the lateral half-vortices of opposite sign are separated by double arrows. The same set 
of rules for extracting the distribution of % from the distribution of half-vortices applies 
also to all the other states with the same energy discussed below. 

The energy of the considered state (calculated per triple site of dice lattice) is given by 



lattice |gjTT|. 



E = V(9 1 )+V(9 2 ) + V(9 3 ), 



(7) 



whereas general constraints (^) are reduced to 



20! + 29 3 = tx 



(8) 



for the central half-vortex of each cluster and 



- 9 X + 29 2 + 9 3 = 7T 



(9) 
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for all the other (lateral) half- vortices. Variation of Eq. (0) with constraints @-@ gives 

V'(9 1 ) + V'(9 2 )=V'(6 3 ), (10) 

which (not unexpectedly) coincides with the condition of the current conservation for each 
of the triple sites. The current conservation at each of the six-link sites in the state of Fig. 
3a is ensured automatically (by symmetry). 

For V(8) of the form (|3|) (corresponding to Josephson junction array) the solution of 
Eqs. (D-Qig) gives 

0! = arctan Z^zl ~ io° 2 = - + 1 » 55°, 3 = - - B x » 80°; (11) 
v /2 + 1 4 2 v ; 

whereas for V(6) oc # 2 (the case of a thin wire network) 

Q x = 15°, 6 2 = 60°, # 3 = 75°. (12) 

Thus the values of 9 a are only weakly sensitive to the type of superconducting system which 
manifests itself in the form of the current-phase relation. 



IV. THE ADDITIONAL DEGENERACY 

The regular state depicted in Fig. 3a (Fig. 2a) allows for the construction of zero energy 
domain wall. Fig. 3c shows how one can rearrange the arrows in the lower half of Fig. 3a 
without invalidating constraints (|5|) or current conservation relations, obtaining in such way 
another extremum with the same energy. The same state is shown in Fig. 2b in terms of 
distribution of half- vortices. It looks like a domain wall separating the state (a) (i. e. the state 
shown in Fig. 2a) from another version of the same state in which the triads of negative 
half- vortices change their orientation by 60°. 

It is possible to construct such domain wall on each horizontal line similar to the line 
shown in Fig. 2b. This increases the degeneracy by the factor 2 N (where N is the number of 
available positions of the domain walls) due to the binary possibility of having or not having 
a domain wall at each available position. The regular state constructed by inserting into the 
state of Fig. 2a a domain wall at each available position is shown in Fig. 2c. This state is 
also periodic, but has the higher (24-fold) degeneracy than the state (a). 

Insofar we have discussed only the zero-energy domain walls which are parallel to the 
triads of positive half-vortices. It follows from the symmetry considerations that analogous 
domain walls can be also constructed in parallel to the triads of negative half- vortices. 
However, the energy remains the same only if all domain walls of this type have the same 
orientation, therefore the total increase of the degeneracy due to a possible creation of zero- 
energy domain walls of the type (b) is given by a factor 2 N+1 . Analogous restriction for the 
creation of zero-energy domain walls appears in the case of the frustrated W-model with 
a triangular lattice and / = l/4or/ = l/3 |jl2| . 

In Fig. 3 the central half-vortex of each triad is marked by square brackets. It is not hard 
to notice that all the other (lateral) half-vortices form the rows of alternating pluses and 
minuses. These rows are straight for the regular state of Fig. 3a (Fig. 2a), but the presence 
of domain walls of the type (b) makes them bend (in parallel to each other). 



5 



It turns out possible to interchange pluses and minuses in any of these rows by interchang- 
ing the single and triple arrows on the links which separate the row from the neighboring 
central half- vortices and reversing the double arrows on all the links inside the row. This pro- 
cedure does not invalidate current conservation at any site and does not change the energy 
of the system. 

Such sign reversal in the rows of alternating lateral half-vortices allows to construct the 
zero-energy domain wall of different type shown in Fig. 2d which [like the domain wall of the 
type (b)] also separates two different versions of the state (a), but now with interchanged 
orientations of positive and negative triads. A regular repetition of such domain wall at each 
available position leads to the periodic state shown in Fig. 2e, which like the state (a) has 
the 12-fold degeneracy. 

The zero-energy domain walls of different types can cross each other (as is shown in Fig. 
2f) without increasing the energy of the system. Each time time the wall of type (d) crosses 
the wall of the type (b) it has to change its orientation by 60°, all the walls of the type (d) 
being parallel to each other in each strip between the walls of the type (b). 

Thus the system simultaneously with an arbitrary number of the walls of the type (b) can 
contain also an arbitrary number of the walls of type (d). The total additional degeneracy 
[in comparison with our reference state (a)] is given by the factor 2 N+AI+1 , where M is 
the number of positions available for domain walls of the type (d), that is the number of 
rows of the alternating lateral half- vortices. The most dense network of zero-energy domain 
walls of both types produces the periodic state shown in Fig. 2g, which like the state (c) is 
characterized by the 24-fold degeneracy. 



V. DISCUSSION 

The zero-energy domain walls of the two types described above produce the contribu- 
tion to residual (zero-temperature) entropy which grows with the increase of the system as 
(N + M) In 2, that is slower than its area (which is proportional to NM). I. e. they do not 
lead to appearance of the extensive residual entropy in contrast to the case of the antifer- 
romagnetic Xl^-model with a Kagome lattice in which the manifold of the ground states 
allows for construction of zero-energy domain walls which can form independent closed loops 



of arbitrary length [^,|T^]. This property of the antiferromagnetic XY- model with Kagome 



lattice leads both to a finite extensive entropy [13,14| and to the presence of a hierarchical 



sequence of barriers which may explain the experimentally observed glass-like dynamics 



of the antiferromagnet with such structure |I5|. The family of the states considered in this 
work demonstrates less developed degeneracy (analogous to that encountered in / = 1/3 
XY-model with triangular lattice ||12|| ). 



Since different values of 6 a correspond to different values of V"(6), at finite temperatures 
the accidental degeneracy related to possible formation of zero-energy domain walls can be 
expected to be removed due to the difference in free energy of the small amplitude continuous 
fluctuations (spin waves) in the same way as it happens in the XF-model with triangular 
lattice and / = l/4or/ = l/3 jl2| . Most probably the spin wave contribution to free energy 
will be minimal for the one of the periodic states shown in Fig. 2, which therefore will be 
dominant in the low temperature limit. 
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The zero energy domain walls separating from each other the different versions of this 
state will then acquire a positive free energy, so the phase transition associated with their 
proliferation (and vortex pattern disordering) can be expected to happen only at finite 
temperatures. However in the thin wire networks with almost harmonic interaction function 
V(8) the effects related with the differences in spin wave free energy will be extremely weak 
and therefore the disordering of vortex pattern due to proliferation of domain walls may 
happen already at rather low temperatures. 

The strong disordering of vortex pattern observed in FF superconducting network with 
dice lattice geometry J/J can also have some relation to geometrical irregularities. Gupta and 
Teitel have recently shown that in the case of the FF XY-model with square lattice the 
irregularities of so-called "positional disorder" type (uncorrelated lattice sites displacements, 
etc.) produce an effective random field for the Ising-type variables M (describing the signs of 
the half-vortices) and therefore induce the destruction of long-range order (at large enough 
scales) even if the disorder is small. In the system which allows for formation of the zero- 
energy domain walls the relevance of this mechanism may be strongly amplified. 
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FIGURES 






FIG. 1. Dice lattice is periodic and has hexagonal symmetry. It consists of the sites with coor- 
dination numbers 3 and 6. All elementary plaquettes are rhombic. 
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g) 

FIG. 2. Filled (empty) circles designate positive (negative) half- vortices. Structures a), c), e) 
and g) are periodic, whereas b), d) and f) include zero-energy domain walls separating different 
periodic states. All states shown have the same energy. 
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FIG. 3. a) Phase representation of the periodic state shown in Fig. 2a; b) elementary pattern, 
repetition of which allows to construct this state; c) phase representation of a zero-energy domain 
wall (shown in Fig. 2b). Three types of arrows correspond to three different values of 
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